Method for initial alignment of radar assisted airborne strapdown inertial navigation system

ABSTRACT

The invention provides a method for initial alignment of radar assisted airborne strapdown inertial navigation system. By calculating the slant distance and angular position between the radar and the airborne inertial navigation equipment, a nonlinear measurement equation for the initial alignment of the radar assisted inertial navigation system is obtained. The unscented Kalman filter algorithm is used to estimate and compensate the error amount of strapdown inertial navigation system to complete the initial alignment task. The significance of the present invention is to provide an in-flight initial alignment solution when the global positioning system is limited, which has fast convergence speed and high estimation accuracy and has high engineering application value.

TECHNICAL FIELD

The invention relates to an alignment method, in particular to a method for initial alignment of radar assisted airborne strapdown inertial navigation system.

BACKGROUND TECHNOLOGY

The initial alignment technology of strapdown inertial navigation system is the key technology of inertial navigation, which directly affects the accuracy of inertial navigation. The current research on aerial alignment is mostly focused on the alignment of shipborne weapons, but the airborne strapdown inertial navigation system also needs to be restarted in special operations. At present, the commonly used in-flight alignment method is Global Positioning System (GPS) assisted alignment. Considering that GPS is susceptible to blockade restrictions and unavailability in wartime, it is particularly important to find other auxiliary methods to achieve SINS initial alignment under special circumstances. Its ranging range can reach thousands to tens of thousands of kilometers, the precision of angle measurement is high, and it has the characteristics of continuous tracking, high precision measurement and high data rate output. The setting is simple, which can be provided to the onboard navigation system via wireless transmission. At present, there is not much research on the initial alignment of radar-assisted inertial navigation equipment. Usually, the position information in the earth coordinate system output by the radar is directly used to construct the quantity measurement. In fact, this model is inaccurate. As the distance between the target and the radar increases, its linearized position error will become larger, that is, the measurement noise will also change. This is very unfavorable for the state estimation and usually causes a large alignment error.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide a method for initial alignment of radar assisted airborne strapdown inertial navigation system, which is a method of initial alignment in airborne strapdown inertial navigation system assisted with slant distance and angular position information measurement by radar.

The purpose of the present invention is achieved as follows:

Step 1: position and track the target aircraft through the tracking radar configured on the ship;

Step 2: tracking radar measures the slant distance, angular position of the target aircraft and the position of the radar will be provided to the airborne inertial navigation system through wireless transmission;

Step 3: construct the initial alignment state equation of strapdown inertial navigation system;

Step 4: construct a nonlinear measurement equation for strapdown inertial navigation system initial alignment;

Step 5: use the unscented Kalman filter to estimate and compensate the inertial navigation error.

The invention also includes such structural features:

1. Step 3 specifically includes:

A state X is selected as:

X=[ϕ^(T) (δv ^(n))^(T)(δP)^(T)(ε^(b))^(T)(∇^(b))^(T)]^(T)

Among them: n is the navigation coordinate system, which coincides with the local geographic coordinate system, and its x, y, and z axes point to east, north, and vertical, respectively; b is the carrier coordinate system, and its x, y, and z point to the right, front and top of the carrier, respectively; δP=[δL δλ≢h]^(T) is the position error vector, δL is the latitude position error, δλ is the longitude position error, δh height position error; δv^(n)=[δv_(E) δv_(N) δv_(v)]^(T) is the speed error vector, δv_(E) is the eastward velocity error, δv_(N) is the northward velocity error, δv_(u) is the celestial velocity error; ϕ=[ϕ_(e) ϕ_(n) ϕ_(u)]^(T) is the platform misalignment angle vector, φ_(e), Φ_(n) and Φ_(u) are the platform misalignment angles in the east, north, and vertical directions, respectively, ε^(b)=[ε_(x) ε_(y) ε_(z)]^(T) is the constant drift vector of the gyro, ε_(x), ε_(y) and ε_(z) are the constant drifts of the gyro in the x, y, and z axes, respectively, ∇^(b)=[∇_(x) ∇_(y) ∇_(z)]^(T) is the accelerometer constant bias vector, and ∇_(x), ∇_(y) and ∇_(z) are the accelerometer constant bias in the x, y, and z axes;

According to the selected state parameters, the initial alignment state equation of the inertial navigation system is

$\quad\left\{ \begin{matrix} {\overset{.}{\varphi} = {{{- \omega_{in}^{n}} \times \varphi} + {M_{12}\delta \; v^{n}} + {M_{13}\delta \; P} - {C_{b}^{n}\left( {ɛ^{b} + ɛ_{w}^{b}} \right)}}} \\ {{\delta \; {\overset{.}{v}}^{n}} = {{\left( {C_{b}^{n}f^{b}} \right) \times \varphi} + {\left( {{v^{n} \times M_{12}} - \left\lbrack {\left( {{2\omega_{ie}^{n}} + \omega_{en}^{n}} \right) \times} \right\rbrack} \right)\delta \; v^{n}} + {M_{23}\delta \; P} + {C_{b}^{n}\left( {\nabla^{b}{+ \nabla_{w}^{b}}} \right)}}} \\ {{\delta \; \overset{.}{P}} = {{M_{32}\delta \; v^{n}} + {M_{33}\delta \; P}}} \\ {{\overset{.}{ɛ}}^{b} = 0} \\ {{\overset{.}{\nabla}}^{b}{= 0}} \end{matrix} \right.$

Among them: ω_(in) ^(n)=ω_(ie) ^(n)+ω_(en) ^(n) is the projection of rotation angular velocity of the navigation system relative to the inertial system in the navigation system, which includes two vector parts: ω_(ie) ^(n) is the projection vector of the earth's rotation angular velocity in the navigation coordinate system, ω_(en) ^(n) is the projection vector of the rotation angular velocity of the navigation system relative to the inertial system caused by the motion of the carrier on the surface of the earth in the navigation coordinate system.

${\omega_{ie}^{n} = \begin{bmatrix} 0 & {\omega_{ie}\cos \; L} & {\omega_{ie}\sin \; L} \end{bmatrix}^{T}},{\omega_{en}^{n} = \left\lbrack {\frac{- v_{N}}{R_{h}}\ \frac{v_{E}}{R_{h}}\ \frac{v_{E}\tan \; L}{R_{h}}} \right\rbrack^{T}}$

and ω_(ie) the earth's rotation angular rate scalar, L is the local latitude, R_(h) is the distance between the carrier and the center of the earth, where R_(h)=R_(e)+h, R_(e) is the radius of the earth, h is the altitude of the carrier , v^(n)=[v_(E) v_(N) v_(u)]^(T) is the projection of the carrier velocity vector in the navigation coordinate system. v_(E), v_(N), and v_(U) are the east velocity, north velocity, and vertical velocity, respectively. f^(b) is the specific force vector of the accelerometer output carrier coordinate system,

$\mspace{20mu} {{M_{12} = \begin{bmatrix} 0 & {{- 1}/R_{h}} & 0 \\ {1/R_{h}} & 0 & 0 \\ {\tan \; {L/R_{h}}} & 0 & 0 \end{bmatrix}},\mspace{20mu} {M_{13} = {\begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\sin \; L} & 0 & 0 \\ {\omega_{ie}\cos \; L} & 0 & 0 \end{bmatrix} + \begin{bmatrix} 0 & 0 & {v_{N}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\; \tan \; {L/R_{h}^{2}}} \end{bmatrix}}},{M_{23} = {\left( {v^{n} \times} \right)\left( {{2 \cdot \begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\; \sin \; L} & 0 & 0 \\ {\omega_{ie}\; \cos \; L} & 0 & 0 \end{bmatrix}} + \begin{bmatrix} 0 & 0 & {{- v_{N}}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\; \tan \; {L/R_{h}^{2}}} \end{bmatrix}} \right)}},\mspace{20mu} {M_{32} = \begin{bmatrix} 0 & {1/R_{h}} & 0 \\ {\sec \; {L/R_{h}}} & 0 & 0 \\ 0 & 0 & 1 \end{bmatrix}},\mspace{20mu} {M_{33} = {\begin{bmatrix} 0 & 0 & {{- v_{N}}/R_{h}^{2}} \\ {v_{E}\; \sec \; L\; \tan \; {L/R_{h}}} & 0 & {{- v_{E}}\; \sec {L/R_{h}^{2}}} \\ 0 & 0 & 0 \end{bmatrix}\mspace{14mu} {and}}}}\mspace{14mu}$ $\mspace{20mu} {C_{b}^{n} = \begin{bmatrix} T_{11} & T_{12} & T_{13} \\ T_{21} & T_{22} & T_{23} \\ T_{31} & T_{32} & T_{33} \end{bmatrix}}$

are the attitude matrix of the sub-inertial navigation, T₁₁, T₁₂, T₁₃, T₂₁, T₂₂, T₂₃, T₃₁, T₃₂ and T₃₃ are the elements of the attitude matrix, ε_(w) ^(b)=[ε_(wx) ^(b) ε_(wy) ^(b) ε_(wz) ^(b)]^(T) the Gaussian white noise vector measured by the gyro, ε_(wx) ^(b), ε_(wy) ^(b) and ε_(wz) ^(b) are the x, y and z axis gyro measurement Gaussian white noise, ∇_(w) ^(b)=[∇_(wx) ^(b) ∇_(wy) ^(b) ∇_(wz) ^(b)]^(T) is the white Gaussian vector of accelerometer measurement, and ∇_(wx) ^(b), ∇_(wy) ^(b) and ∇_(wz) ^(b) are the white Gaussian of x, y and z axial acceleration measurement.

2. Step 4 includes:

Measurement Z=[R β α]^(T) includes slant distance R, azimuth angle β and pitch angle α;

Among them:

${Z = {\begin{bmatrix} R \\ \beta \\ \alpha \end{bmatrix} = \begin{bmatrix} \sqrt{\left( {dx^{n}} \right)^{2} + \left( {dy^{n}} \right)^{2} + \left( {dz}^{n} \right)^{2}} \\ {\arctan \frac{{dz}^{n}}{\sqrt{\left( {dx}^{n} \right)^{2} + \left( {dy}^{n} \right)^{2}}}} \\ {{arc}\; \tan \frac{{dx}^{n}}{{dy}^{n}}} \end{bmatrix}}},{\begin{bmatrix} {dx}^{n} & {dy}^{n} & {dz}^{n} \end{bmatrix}^{T} = {{C_{p_{o}}^{n}\begin{bmatrix} {dx}^{e} & {dy}^{e} & {dz}^{e} \end{bmatrix}}^{T}\mspace{14mu} {{and}\begin{bmatrix} {dx}^{n} & {dy}^{n} & {dz}^{n} \end{bmatrix}}T}}$

are the relative position vector of the target and the radar in the navigation coordinate system, C_(P) _(o) ^(n) is the coordinate conversion matrix between the earth's rectangular coordinate system and the navigation coordinate system. The coordinate of the Earth's rectangular coordinate system is

${\begin{bmatrix} {dx^{e}} \\ {dy^{e}} \\ {dz^{e}} \end{bmatrix} = {\begin{bmatrix} {\left( {R_{e} + h_{p}} \right){\cos \left( L_{p} \right)}{\cos \left( \lambda_{p} \right)}} \\ {\left( {R_{e} + h_{p}} \right){\cos \left( L_{p} \right)}{\sin \left( \lambda_{p} \right)}} \\ {\left( {R_{e} + h_{p}} \right){\sin \left( L_{p} \right)}} \end{bmatrix} - \begin{bmatrix} x_{p_{o}}^{e} \\ y_{p_{o}}^{e} \\ z_{p_{o}}^{e} \end{bmatrix}}},$

e represents the Earth's rectangular coordinate system. L_(p)=L_(p) ^(s)−δL is the true latitude, λ_(p)=λ_(p) ^(s)−δλ is true longitude, h_(p)=h_(p) ^(s)−δh is true altitude, and L_(p) ^(s), λ_(p) ^(s) and h_(p) ^(s) are the position resolved by the inertial navigation system.

Then the measurement equation for the initial alignment of the radar-assisted strapdown inertial navigation system is:

$Z = {{H\left( {{\delta L},{\delta \lambda},{\delta h}} \right)} + \begin{bmatrix} \omega_{R} \\ \omega_{\beta} \\ \omega_{\alpha} \end{bmatrix}}$

Among them: ω_(R), ω_(α) and ω_(β) top are white noises that conform to the zero-mean Gaussian distribution, and the expression of the nonlinear function H can be obtained by the above substitution.

Compared with the prior art, the beneficial effects of the present invention are: the present invention is based on the slant distance and angular position provided by the radar, considering the transfer relationship between the positioning error of the strapdown inertial navigation system and the slant distance and angular position, the present invention proposes a new alignment scheme for the measurement model, which uses the slant distance and angular position as measurement information to achieve alignment. First, the present invention provides a new solution for the initial alignment of the airborne inertial navigation system when the global positioning system is blocked, and has high engineering application value. Second, compared with the traditional radar-assisted inertial navigation system initial alignment scheme, the advantages of the present invention are reflected in the fact that the position coordinates after the radar measurement parameters are linearized are not selected as the measurement, which avoids the problem of the statistical characteristics of the measurement noise changing with the distance. The invention directly uses the slant distance and angular position information obtained by radar measurement as the quantity measurement, and makes full use of the original measurement information. Since the statistical characteristics of the measurement noise meet the requirements of the optimal estimation, the Kalman filter method can be used to estimate the optimal status. At the same time, compared with the existing alignment scheme, the proposed scheme can complete the high-precision initial alignment task in a larger distance range.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1 is a schematic diagram of the present invention;

FIG. 2 is a description of the radar target measurement parameters of the present invention;

FIG. 3a is a comparison diagram of horizontal attitude error between the present invention and the existing scheme under the condition 1 in simulation experiment, FIG. 3b is a comparison diagram of azimuth attitude error between the present invention and the existing scheme under the condition 1 in simulation experiment.

FIG. 4a is a comparison diagram of horizontal attitude error between the present invention and the existing scheme under the condition 2 in simulation experiment, FIG. 4b is a comparison diagram of azimuth attitude error between the present invention and the existing scheme under the condition 2 in simulation experiment.

DETAILED DESCRIPTION OF THE INVENTION

The present invention will be further described in detail below with reference to the drawings and specific embodiments.

As shown in FIG. 1, it is a schematic diagram of an initial alignment solution for radar-assisted strapdown inertial navigation system provided by the present invention. The measurement noise statistical characteristics vary with the slant distance. The invention has the advantages of high alignment accuracy and short alignment time. It includes the following steps:

Step 1: the tracking radar is configured on the ship to locate and track the target aircraft. As shown in FIG. 2, P_(o) is the location of the radar carrier, and the polar coordinates P (R, α, β) of any target P in the air are measured with it as the origin. R is called the slant distance, which is the linear distance from the radar to the target. α represents the azimuth angle, α represents the azimuth angle, which is the angle between PoB projected on the horizontal plane between the radar and the target line PoP and the true north direction on the horizontal plane. β represents the elevation angle, which is the angle between the PoB projected on the horizontal plane between the radar and the target line PoP and the lead vertical plane. It is also called the inclination angle or the height angle.

Step 2: tracking radar measures the slant distance, angular position of the target aircraft and the location of the radar carrier will be provided to the airborne inertial navigation system through wireless transmission. The location P_(o) of the radar carrier is provided by a high-precision inertial navigation device, and its positioning accuracy is high, so the influence of its position error can be ignored. Because of its high positioning accuracy, its position error can be ignored.

Step 3: construct the initial alignment state parameters and state equations of the strapdown inertial navigation system. The state quantity is selected as:

X=[ϕ^(T) (δv ^(n))^(T)(δP)^(T)(∇^(b))^(T)]^(T)

n is the navigation coordinate system, which coincides with the local geographic coordinate system, and its x, y, and z axes point to east, north, and vertical, respectively; b is the carrier coordinate system, and its x, y, and z point to the right , front and top of carrier; δP=[δL δλ δh]^(T) is the position error vector, δL is the latitude position error, δλ is the longitude position error, δh height position error; δv^(n)=[δv_(E) δv_(N) δv_(U)]^(T)is the speed error vector , δv_(E) is the eastward velocity error, δv_(N) is the northward velocity error, δv_(U) is the celestial velocity error; ϕ=[ϕ_(e) ϕ_(n) ϕ_(u)]^(T) is the platform misalignment angle vector, Φ_(e), Φ_(n) and Φ_(u) are the platform misalignment angles in the east, north, and vertical directions, respectively, ε^(b)=[ε_(x) ε_(y) ε_(z)]^(T) is the constant drift vector of the gyro, ε_(x), ε_(y) and ε_(z) are the constant drifts of the gyro in the x, y, and z axes, respectively, ∇^(b)=[∇_(x) ∇_(y) ∇_(Z)]^(T) is the accelerometer constant bias vector, and ∇_(x), ∇_(y),∇_(z) is the accelerometer constant bias in the x, y, and z axes, and T represents transpose.

Further, according to the selected state parameters can be obtained inertial navigation system initial alignment state equation is:

$\quad\left\{ \begin{matrix} {\overset{.}{\varphi} = {{{- \omega_{in}^{n}} \times \varphi} + {M_{12}\delta \; v^{n}} + {M_{13}\delta \; P} - {C_{b}^{n}\left( {ɛ^{b} + ɛ_{w}^{b}} \right)}}} \\ {{\delta \; {\overset{.}{v}}^{n}} = {{\left( {C_{b}^{n}f^{b}} \right) \times \varphi} + {\left( {{v^{n} \times M_{12}} - \left\lbrack {\left( {{2\omega_{ie}^{n}} + \omega_{em}^{n}} \right) \times} \right\rbrack} \right)\delta \; v^{n}} + {M_{23}\delta \; P} + {C_{b}^{n}\left( {\nabla^{b}{+ \nabla_{w}^{b}}} \right)}}} \\ {{\delta \; \overset{.}{P}} = {{M_{32}\delta \; v^{n}} + {M_{33}\delta \; P}}} \\ {{\overset{.}{ɛ}}^{b} = 0} \\ {{\overset{.}{\nabla}}^{b}{= 0}} \end{matrix} \right.$

Among them, the point on the state quantity represents the first derivative .ω_(in) ^(n)=ω_(ie) ^(n)+ω_(en) ^(n) is the projection of rotation angular velocity of the navigation system relative to the inertial system in the navigation system, which includes two vector parts: ω_(ie) ^(n) is the projection vector of the earth's rotation angular velocity in the navigation coordinate system, ω_(en) ^(n) is the projection vector of the rotation angular velocity of the navigation system relative to the inertial system caused by the motion of the carrier on the surface of the earth in the navigation coordinate system.

${\omega_{ie}^{n} = \begin{bmatrix} 0 & {\omega_{ie}\cos \; L} & {\omega_{ie}\sin \; L} \end{bmatrix}^{T}},{\omega_{en}^{n} = \left\lbrack {\frac{- v_{N}}{R_{h}}\ \frac{v_{E}}{R_{h}}\ \frac{v_{E}\tan \; L}{R_{h}}} \right\rbrack^{T}}$

and ω_(ie) are the earth's rotation angular rate scalar, L is the local latitude, R_(h) is the distance between the carrier and the center of the earth, where R_(h)=R_(e)+h, R_(e) is the radius of the earth, h is the altitude of the carrier, v_(n)=[v_(E) v_(N) v_(U)]^(T) is the projection of the carrier velocity vector in the navigation coordinate system. v_(E), v_(N), and v_(U) are the east velocity, north velocity, and vertical velocity, respectively. f^(b) is the specific force vector of the accelerometer output carrier coordinate system,

$\mspace{20mu} {{M_{12} = \begin{bmatrix} 0 & {{- 1}/R_{h}} & 0 \\ {1/R_{h}} & 0 & 0 \\ {\tan \; {L/R_{h}}} & 0 & 0 \end{bmatrix}},\mspace{20mu} {M_{13} = {\begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\; \sin \; L} & 0 & 0 \\ {\omega_{ie}\cos \; L} & 0 & 0 \end{bmatrix} + \begin{bmatrix} 0 & 0 & {v_{N}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\tan {L/R_{h}^{2}}} \end{bmatrix}}},{M_{23} = {\left( {v^{n} \times} \right)\left( {{2 \cdot \begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\sin \; L} & 0 & 0 \\ {\omega_{ie}\cos \; L} & 0 & 0 \end{bmatrix}} + \ \begin{bmatrix} 0 & 0 & {v_{N}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\tan \; {L/R_{h}^{2}}} \end{bmatrix}} \right)}},\mspace{20mu} {M_{32} = \begin{bmatrix} 0 & {1/R_{h}} & 0 \\ {\sec \; {L/R_{h}}} & 0 & 0 \\ 0 & 0 & 1 \end{bmatrix}},\mspace{20mu} {M_{33} = {\begin{bmatrix} 0 & 0 & {{- v_{N}}/R_{h}^{2}} \\ {v_{E}\sec \; L\; \tan \; {L/R_{h}}} & 0 & {{- v_{E}}\sec \; {L/R_{h}^{2}}} \\ 0 & 0 & 0 \end{bmatrix}\mspace{14mu} {and}}}}\mspace{14mu}$ $\mspace{20mu} {C_{b}^{n} = \begin{bmatrix} T_{11} & T_{12} & T_{13} \\ T_{21} & T_{22} & T_{23} \\ T_{31} & T_{32} & T_{33} \end{bmatrix}}$

are the attitude matrix of the sub-inertial navigation, T₁₁, T₁₂, T₁₃, T₂₁, T₂₂, T₂₃, T₃₁, T₃₂ and T₃₃ are the elements of the attitude matrix, ε_(w) ^(b)=[ε_(wx) ^(b) ε_(wy) ^(b) ε_(wz) ^(b)]^(T) is the Gaussian white noise vector measured by the gyro, ε_(wx) ^(b), ε_(wy) ^(b) and ε_(wz) ^(b) are the x, y and z axis gyro measurement Gaussian white noise, ∇_(w) ^(b)=[∇_(wx) ^(b) ∇_(wy) ^(b) ∇_(wz) ^(b)]^(T) is the white Gaussian vector of accelerometer measurement, and ∇_(wx) ^(b), ∇_(wy) ^(b) and ∇_(wz) ^(b) are the white Gaussian of x, y and z axial acceleration measurement.

Step 4: select the initial alignment measurement parameters of the strapdown inertial navigation system and construct a nonlinear measurement equation. In the existing initial alignment method of the radar-assisted strapdown inertial navigation system, the position error is selected as the measurement. In the present invention, the slant distance and the angular position are used as the measurement; then the measurement: Z=[R β 60 ]^(T) is respectively slant distance R, azimuth angle β and pitch angle α.

Among them:

${Z = {\begin{bmatrix} R \\ \beta \\ \alpha \end{bmatrix} = \begin{bmatrix} \sqrt{\left( {dx^{n}} \right)^{2} + \left( {dy^{n}} \right)^{2} + \left( {dz}^{n} \right)^{2}} \\ {\arctan \frac{{dz}^{n}}{\sqrt{\left( {dx}^{n} \right)^{2} + \left( {dy}^{n} \right)^{2}}}} \\ {{arc}\; \tan \frac{{dx}^{n}}{{dy}^{n}}} \end{bmatrix}}},{\begin{bmatrix} {dx}^{n} & {dy}^{n} & {dz}^{n} \end{bmatrix}^{T} = {C_{p_{o}}^{n}\begin{bmatrix} {dx}^{e} & {dy}^{e} & {dz}^{e} \end{bmatrix}}^{T}}$

and [dx^(n) dy^(n) dz^(n)]^(T) are the relative position vector of the target and the radar in the navigation coordinate system, C_(p) _(o) ^(n) is the coordinate conversion matrix between the earth's rectangular coordinate system and the navigation coordinate system. The coordinate of the earth's rectangular coordinate system is

${\begin{bmatrix} {dx^{e}} \\ {dy^{e}} \\ {dz^{e}} \end{bmatrix} = {\begin{bmatrix} {\left( {R_{e} + h_{p}} \right)\cos \; \left( L_{p} \right)\cos \; \left( \lambda_{p} \right)} \\ {\left( {R_{e} + h_{p}} \right)\cos \; \left( L_{p} \right)\sin \; \left( \lambda_{p} \right)} \\ {\left( {R_{e} + h_{p}} \right)\sin \; \left( L_{p} \right)} \end{bmatrix} - \begin{bmatrix} x_{p_{o}}^{e} \\ y_{p_{o}}^{e} \\ z_{p_{o}}^{e} \end{bmatrix}}},$

e represents the earth's rectangular coordinate system. L_(p)=L_(p) ^(s)−δL is the true latitude, λ_(p)=λ_(p) ^(s)−δλ0 is true longitude, h_(p)=h_(p) ^(s)−δh is true altitude, and L_(p) ^(s), λ_(p) ^(s) and h_(p) ^(s) are the position resolved by the inertial navigation system. Then the measurement equation for the initial alignment of the radar-assisted strapdown inertial navigation system is:

$Z = {{H\left( {{\delta L},{\delta\lambda},{\delta \; h}} \right)} + {\begin{bmatrix} \omega_{R} \\ \omega_{\beta} \\ \omega_{\alpha} \end{bmatrix}.}}$

${Z = {{H\left( {{\delta L},{\delta \lambda},{\delta h}} \right)} + \begin{bmatrix} \omega_{R} \\ \omega_{\beta} \\ \omega_{\alpha} \end{bmatrix}}},$

ω_(R), ω_(α) and ω_(β) are white noises that conform to a zero-mean Gaussian distribution, and the expression of the nonlinear function H can be obtained by the above substitution.

Step 5: use the unscented Kalman filter to estimate and compensate the strapdown inertial navigation system error.

The system equations and measurement equations for the initial alignment of the radar-assisted strapdown inertial navigation system are given in steps 3 and 4. The initial alignment task can be completed only by estimating and compensating the state quantities. Because the measurement equation is non-linear, this scheme uses the unscented Kalman filter algorithm for state estimation.

(1) Select the initial filter value

{circumflex over (X)}₀=EX₀

P _(o) =E[(X ₀-{circumflex over (X)}₀)][(X ₀-{circumflex over (X)}₀)^(T)]

System dimension n=15

The weights are:

${W_{0}^{(m)} = \frac{\lambda}{n + \lambda}},\mspace{11mu} {W_{0}^{(c)} = {\frac{\lambda}{n + \lambda} + 1 - a^{2} + b}}\;,\mspace{14mu} {W_{i}^{(m)} = {W_{i}^{(c)} = \frac{\lambda}{2\left( {n + \lambda} \right)}}}$      i = 1, 2, …  , 2n

γ=√{square root over (n+λ)}, λ=α²(n+κ)−n a is a very small positive number, 10⁻⁴≤α≤1, κ=3-n, b=2 can be selected.

(2) Calculate 2n+1 σ samples when k-1 (k=1, 2, 3, . . . )

{tilde over (X)}_(k-1) ⁽⁰⁾={circumflex over (X)}_(k-1)

{tilde over (X)}_(k-1) ^((i))={circumflex over (X)}_(k-1)+γ(√{square root over (P _(k-1))})_((i)) i=1,2, . . . n

{tilde over (X)}_(k-1) ^((i))={circumflex over (X)}_(k-1)−γ(√{square root over (P _(k-1))})_((i-n)) i=n+1, n+2, . . . 2n

(3) A predictive model for computing k time

${\chi_{{k/k} - 1}^{*{(i)}} = {{{f\left\lbrack {{\overset{\sim}{\chi}}_{k - 1}^{(i)},u_{k - 1}} \right\rbrack}\mspace{31mu} i} = 0}},1,2,\ldots \mspace{14mu},{2n}$ ${\hat{X}}_{{k/k} - 1} = {\sum\limits_{i = 0}^{2n}{W_{i}^{(m)}\chi_{{k/k} - 1}^{*{(i)}}}}$ $P_{{k/k} - 1} = {{\sum\limits_{i = 0}^{2n}{{W_{i}^{(c)}\left\lbrack {\chi_{{k/k} - 1}^{*{(i)}} - {\hat{X}}_{{k/k} - 1}} \right\rbrack}\;\left\lbrack {\chi_{{k/k} - 1}^{*{(i)}} - {\hat{X}}_{{k/k} - 1}} \right\rbrack}^{T}} + Q_{k - 1}}$

(4) Calculate the one-step prediction sample point at the time of k)

X _(k/k-1) ⁽⁰⁾={circumflex over (X)}_(k/k-1)

X _(k/k-1) ^((i))={circumflex over (X)}_(k/k-1)+γ(√{square root over (P _(k/k-1))})_((i)) i=1,2, . . . , n

X _(k/k-1) ^((i))={circumflex over (X)}_(k/k-1)−γ(√{square root over (P _(k/k-1))})_((i-n)) i=n+1, n+2, . . . , 2n

(5) Calculation P_((XZ)k/k-1), P_((ZZ)k/k-1)

Z_(k/k − 1)^((i)) = h[χ_(k/k − 1)^((i))],  i = 0, 1, 2, …  , 2n ${\overset{\hat{}}{Z}}_{{k/k} - 1} = {{\sum\limits_{i = 0}^{2n}{W_{i}^{(m)}Z_{{k/k} - 1}^{(i)}P_{{({XZ})}_{{k/k} - 1}}}} = {\sum\limits_{i = 0}^{2n}{{W_{i}^{(c)}\left\lbrack {\chi_{{k/k} - 1}^{(i)} - {\overset{\hat{}}{X}}_{{k/k} - 1}} \right\rbrack}\left\lbrack {Z_{{k/k} - 1}^{(i)} - {\overset{\hat{}}{Z}}_{{k/k} - 1}} \right\rbrack}^{T}}}$ $P_{{({ZZ})}_{{k/k} - 1}} = {{\sum\limits_{i = 0}^{2n}{{W_{i}^{(c)}\left\lbrack {Z_{{k/k} - 1}^{(i)} - {\overset{\hat{}}{Z}}_{{k/k} - 1}} \right\rbrack}\left\lbrack {Z_{{k/k} - 1}^{(i)} - {\overset{\hat{}}{Z}}_{{k/k} - 1}} \right\rbrack}^{T}} + R_{k}}$

(6) Calculate the gain matrix

K_(k) = P_((XZ)_(k/k − 1))P_((ZZ)_(k/k − 1))⁻¹

(7) Calculate the filter value

${{\overset{\hat{}}{X}}_{k} = {{\overset{\hat{}}{X}}_{{k/k} - 1} + {K_{k}\left\lbrack {Z_{k} - {\overset{\hat{}}{Z}}_{{k/k} - 1}} \right\rbrack}}}{P_{k} = {P_{{k/k} - 1} - {K_{k}P_{{({ZZ})}_{{k/k} - 1}}K_{k}^{T}}}}$

(8) Through the above process, the navigation error of the strapdown inertial navigation system can be estimated, so as to perform closed-loop correction and complete the initial alignment.

The technical solution of the present invention is simulated and verified in combination with specific values below:

Simulation conditions: the initial position error of the inertial navigation device is set to 10√{square root over (3m)} the parameters of the inertial measurement unit as follows: the constant drift of the gyro is 0.01° h, The random drift is 0.001°/√{square root over (h)}, the accelerometer constant bias is 3×10⁻⁴ g, and the random drift is 5×10⁻⁵ g√{square root over (s)}, the sampling time interval is 10 ms; the slant distance error measured by radar is 10m (1σ), the pitch angle error is 0.1° (1σ), and the azimuth angle error is 0.3° (1σ).

Because the ship's position and attitude are provided by the combination of ship inertial navigation equipment and auxiliary equipment, its position error and attitude error can be ignored. To simplify the simulation complexity, in order to simplify the complexity of simulation, considering that the ship is still and the “true north” provided by the ship is error-free. The height of the radar is 5 m, the update period is 1 s, the filter filtering period is ls, the filter is closed-loop corrected, the simulation time is set to 300 s; the flight speed is 80 m/s, and the flight height is 1000 m.

In order to verify the effectiveness of the present invention, the simulation environment is set: the slant distance of simulation condition 1 is less than 10 km, and the slant distance of simulation condition 2 is more than 50 km. Compared with the existing radar-assisted alignment scheme. FIG. 3a and FIG. 3b are the attitude error comparison diagrams of the present invention and the existing scheme for 200 times Monte Carlo simulation results under conditions 1 in simulation experiment. FIG. 4a and FIG. 4b are the attitude error comparison diagrams of the present invention and the existing scheme for 200 times Monte Carlo simulation results under conditions 2 in simulation experiment. Among them, the thin black solid line is the mean curve of the traditional measurement model scheme, and the thin black dashed line is the 3σ curve of the traditional measurement model scheme. Among them, the thick black solid line is the mean curve of the new measurement model scheme, and the thick black dotted line is the 3σ curve of the new measurement model scheme; Φ_(e), Φ_(n) and Φ_(u) are the pitch error angle, the roll error angle, and the heading error angle, respectively.

In summary, the present invention provides a method for initial alignment of radar assisted airborne strapdown inertial navigation system. By calculating the slant distance and angular position between the radar and the airborne inertial navigation equipment, a nonlinear measurement equation for the initial alignment of the radar assisted inertial navigation system is obtained. The unscented Kalman filter algorithm is used to estimate and compensate the error amount of strapdown inertial navigation system to complete the initial alignment task. The significance of the present invention is to provide an in-flight initial alignment solution when the global positioning system is limited, which has fast convergence speed and high estimation accuracy and has high engineering application value. 

1. A method for initial alignment of radar assisted airborne strapdown inertial navigation system is characterized by the following steps: Step 1: position and track the target aircraft through the tracking radar configured on the ship; Step 2: tracking radar measures the slant distance, angular position of the target aircraft and the position of the radar will be provided to the airborne inertial navigation system through wireless transmission; Step 3: construct the initial alignment state equation of strapdown inertial navigation system; Step 4: construct a nonlinear measurement equation for strapdown inertial navigation system initial alignment; Step 5: use the unscented Kalman filter to estimate and compensate the inertial navigation error.
 2. According to claim 1, a method for initial alignment of radar assisted airborne strapdown inertial navigation system is characterized in that step 3 specifically includes: A state X is selected as: X=[ϕ^(T) (δv ^(n))^(T)(δP)^(T)(ε^(b))^(T)(∇^(b))^(T)]^(T) Among them: n is the navigation coordinate system, which coincides with the local geographic coordinate system, and its x, y, and z axes point to east, north, and vertical, respectively; b is the carrier coordinate system, and its x, y, and z point to the right, front and top of carrier; δP=[δL δλ δh] is the position error vector, δL is the latitude position error, δλ is the longitude position error, δh height position error; δv^(n)=[δv_(E) δv_(N) δv_(U)]^(T); is the speed error vector , δv_(E) is the eastward velocity error, δv_(N) is the northward velocity error, δv_(U) is the celestial velocity error; ϕ=[ϕ_(e) ϕ_(n) ϕ_(u)]^(T) is the platform misalignment angle vector, φ_(e), φ_(n) and Φ_(u) are the platform misalignment angles in the east, north, and vertical directions, respectively, ε^(b)=[ε_(x) ε_(y) ε_(z)]^(T) is the constant drift vector of the gyro, ε_(x), ε_(y) and ε_(z) are the constant drifts of the gyro in the x, y, and z axes, respectively, ∇_(b)=[∇_(x) ∇_(y) ∇_(z)]^(T) is the accelerometer constant bias vector, and ∇_(x), ∇_(y), ∇_(z) is the accelerometer constant bias in the x, y, and z axes; According to the selected state parameters, the initial alignment state equation of the inertial navigation system is $\left\{ {\begin{matrix} {\overset{.}{\varphi} = {{{- \omega_{i\; n}^{n}} \times \varphi} + {M_{12}\delta \; v^{n}} + {M_{13}\delta \; P} - {C_{b}^{n}\left( {ɛ^{b} + ɛ_{w}^{b}} \right)}}} \\ \begin{matrix} {{\delta \; {\overset{.}{v}}^{n}} = {{\left( {C_{b}^{n}f^{b}} \right) \times \varphi} + {\left( {{v^{n} \times M_{12}} - \left\lbrack {\left( {{2\; \omega_{ie}^{n}} + \omega_{en}^{n}} \right) \times} \right\rbrack} \right)\delta \; v^{n}} +}} \\ {{M_{23}\delta \; P} + {C_{b}^{n}\left( {\nabla^{b}{+ \nabla_{w}^{b}}} \right)}} \end{matrix} \\ {{\delta \; \overset{.}{P}} = {{M_{32}\delta \; v^{n}} + {M_{33}\delta \; P}}} \\ {{\overset{.}{ɛ}}^{b} = 0} \\ {{\overset{.}{\nabla}}^{b}{= 0}} \end{matrix}\quad} \right.$ Among them: ω_(in) ^(n)=ω_(ie) ^(n)+ω_(en) ^(n) is the projection of the rotation angular velocity of the navigation system relative to the inertial system on the navigation system, which includes two vector parts: ω_(ie) ^(n) is the projection vector of the earth's rotation angular velocity in the navigation coordinate system, ω_(en) ^(n) is the projection vector of the rotation angular velocity of the navigation system relative to the inertial system caused by the movement of the carrier on the surface of the e_(an)art_(d)h_(oo)i_(i)n_(e) the navigation coordinate system, ${\omega_{ie}^{n} = \begin{bmatrix} 0 & {\omega_{ie}\cos \; L} & {\omega_{ie}\sin \; L} \end{bmatrix}^{T}},{\omega_{en}^{n} = \begin{bmatrix} \frac{- v_{N}}{R_{h}} & \frac{v_{E}}{R_{h}} & \frac{v_{E}\tan \; L}{R_{h}} \end{bmatrix}^{T}}$ and ω_(ie) are the earth's rotation angular rate scalar, L is the local latitude, R_(h) is the distance between the carrier and the center of the earth, where R_(h)=R_(e)+h, R_(e) is the radius of the earth, h is the altitude of the carrier , v^(n)=[v_(E) v_(N) v_(U)]^(T) is the projection of the carrier velocity vector in the navigation coordinate system. v_(E), v_(N), and v_(u) are the east velocity, north velocity, and vertical velocity, respectively. f^(b) is the specific force vector of the accelerometer output carrier coordinate system, $\mspace{20mu} {{M_{12} = \begin{bmatrix} 0 & {{- 1}/R_{h}} & 0 \\ {1/R_{h}} & 0 & 0 \\ {\tan \; {L/R_{h}}} & 0 & 0 \end{bmatrix}},\mspace{20mu} {M_{13} = {\begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\sin \; L} & 0 & 0 \\ {\omega_{ie}\cos \; L} & 0 & 0 \end{bmatrix} + \begin{bmatrix} 0 & 0 & {v_{N}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\; \sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\tan \; {L/R_{h}^{2}}} \end{bmatrix}}},{M_{23} = {\left( {v^{n} \times} \right)\left( {{2 \cdot \begin{bmatrix} 0 & 0 & 0 \\ {{- \omega_{ie}}\sin \; L} & 0 & 0 \\ {\omega_{ie}\cos \; L} & 0 & 0 \end{bmatrix}} + \ \begin{bmatrix} 0 & 0 & {v_{N}/R_{h}^{2}} \\ 0 & 0 & {{- v_{E}}/R_{h}^{2}} \\ {v_{E}\sec^{2}{L/R_{h}}} & 0 & {{- v_{E}}\tan \; {L/R_{h}^{2}}} \end{bmatrix}} \right)}},\mspace{20mu} {M_{32} = \begin{bmatrix} 0 & {1/R_{h}} & 0 \\ {\sec \; {L/R_{h}}} & 0 & 0 \\ 0 & 0 & 1 \end{bmatrix}},\mspace{20mu} {M_{33} = {\begin{bmatrix} 0 & 0 & {{- v_{N}}/R_{h}^{2}} \\ {v_{E}\sec \; L\; \tan \; {L/R_{h}}} & 0 & {{- v_{E}}\sec \; {L/R_{h}^{2}}} \\ 0 & 0 & 0 \end{bmatrix}\mspace{14mu} {and}}}}$ $\mspace{20mu} {C_{b}^{n} = \begin{bmatrix} T_{11} & T_{12} & T_{13} \\ T_{21} & T_{22} & T_{23} \\ T_{31} & T_{32} & T_{33} \end{bmatrix}}$ are the attitude matrix of the sub-inertial navigation, T₁₁, T₁₂, T₁₃, T₂₁, T₂₂, T₂₃, T₃₁, T₃₂ and T₃₃ are the elements of the attitude matrix, ε_(w) ^(b)=[εwx^(b) ε_(wy)ε_(wz) ^(b)]^(T) is the Gaussian white noise vector measured by the gyro, ε_(wx) ^(b), ε_(wy) ^(b) and ε_(wz) ^(b) are the x, y, z axis gyro measurement Gaussian white noise, ∇_(w) ^(b)=[∇_(wx) ^(b) ∇_(wy) ^(b) ∇_(wz) ^(b)]^(T) is the white Gaussian vector of accelerometer measurement, and ∇_(wx) ^(b), ∇_(wy) ^(yb) and ∇_(wz) ^(b) are the white Gaussian of x, y and z axial acceleration measurement.
 3. A method for initial alignment of radar assisted airborne strapdown inertial navigation system according to claim 1 or 2, wherein step 4 specifically includes: Measurement Z=[R β 62 ] ^(T) includes slant distance R, azimuth angle β and pitch angle α; Among them: ${Z = {\begin{bmatrix} R \\ \beta \\ \alpha \end{bmatrix} = \begin{bmatrix} \sqrt{\left( {dx}^{n} \right)^{2} + \left( {dy}^{n} \right)^{2} + \left( {dz}^{n} \right)^{2}} \\ {{arc}\; \tan \frac{{dz}^{n}}{\sqrt{\left( {dx}^{n} \right)^{2} + \left( {dy}^{n} \right)^{2}}}} \\ {{arc}\; \tan \frac{{dx}^{n}}{{dy}^{n}}} \end{bmatrix}}},{\begin{bmatrix} {dx}^{n} & {dy}^{n} & {dz}^{n} \end{bmatrix}^{T} = {C_{p_{o}}^{n}\begin{bmatrix} {dx}^{e} & {dy}^{e} & {dz}^{e} \end{bmatrix}}^{T}}$ and [dx^(n) dy^(n) dz^(n)]^(T) are the relative position vector of the target and the radar in the navigation coordinate system, C_(P) _(o) ^(n) is the coordinate conversion matrix between the earth's rectangular coordinate system and the navigation coordinate system. The coordinate of the earth's rectangular coordinate system is ${\begin{bmatrix} {dx^{e}} \\ {dy^{e}} \\ {dz^{e}} \end{bmatrix} = {\begin{bmatrix} {\left( {R_{e} + h_{p}} \right)\cos \; \left( L_{p} \right)\cos \; \left( \lambda_{p} \right)} \\ {\left( {R_{e} + h_{p}} \right)\cos \; \left( L_{p} \right)\sin \; \left( \lambda_{p} \right)} \\ {\left( {R_{e} + h_{p}} \right)\sin \; \left( L_{p} \right)} \end{bmatrix} - \begin{bmatrix} x_{p_{o}}^{e} \\ y_{p_{o}}^{e} \\ z_{p_{o}}^{e} \end{bmatrix}}},$ e represents the Earth's rectangular coordinate system. L_(p)=L_(p) ^(s)-δL is the true latitude, λ_(p)=λ_(p) ^(s)-δλ is true longitude, h_(p)=h_(p) ^(s)-δh is true altitude, and L_(p) ^(s), λ_(p) ^(s) and h_(p) ^(s) are the position resolved by the inertial navigation system. ${Z = {{H\left( {{\delta L},{\delta \lambda},{\delta h}} \right)} + \begin{bmatrix} \omega_{R} \\ \omega_{\beta} \\ \omega_{\alpha} \end{bmatrix}}},$ Then the measurement equation for the initial alignment of the radar-assisted strapdown inertial navigation system is ${Z = {{H\left( {{\delta L},{\delta \lambda},{\delta h}} \right)} + \begin{bmatrix} \omega_{R} \\ \omega_{\beta} \\ \omega_{\alpha} \end{bmatrix}}},$ Among them: ω_(R), ω_(α)αand ω_(β) are white noises that conform to the zero-mean Gaussian distribution, and the expression of the nonlinear function H can be obtained by the above substitution. 